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> . Abstract 
l> ■ 

£T) ■ We discuss the phenomenology and cosmology of a Standard- like Model inspired by string theory, in 

which the gauge fields are localized on D-branes wrapping certain compact cycles on an underlying 
geometry, whose intersection can give rise to chiral fermions. The energy scale associated with 
\ string physics is assumed to be near the Planck mass. To develop our program in the simplest way, 

CN ' we work within the construct of a minimal model with gauge-extended sector U(3)b x Sp(l)i x 

U (l)i R xU(1)l. The resulting U(l) content gauges the baryon number B, the lepton number L, and 
> a third additional abelian charge Ir which acts as the third isospin component of an SU(2)r. All 

^ , mixing angles and gauge couplings are fixed by rotation of the U(l) gauge fields to a basis diagonal 

in hyper charge Y and in an anomaly free linear combination of Ir and B — L. The anomalous 
Z' gauge boson obtains a string scale Stiickelberg mass via a 4D version of the Green-Schwarz 
mechanism. To keep the realization of the Higgs mechanism minimal, we add an extra SU(2) 
singlet complex scalar, which acquires a VEV and gives a TeV-scale mass to the non- anomalous 
gauge boson Z". The model is fully predictive and can be confronted with dijet and dilepton data 
from LHC8 and, eventually, LHC14. We show that Mz" ~ 3 — 4 TeV saturates current limits 
from the CMS and ATLAS collaborations. We also show that for M z » < 5 TeV, LHC14 will 
reach discovery sensitivity > 5a. After that, we demostrate in all generality that Z" milli-weak 
interactions could play an important role in observational cosmology. Finally, we examine some 
phenomenological aspects of the supersymmetric extension of the D-brane construct. 
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I. INTRODUCTION 



With the turn on of the Large Hadron Collider (LHC) at CERN, a new era of discovery 
has just begun [1-4]. The SU(3) C x SU(2) L x U(1) Y Standard Model (SM) of electroweak 
and strong interactions was once again severely tested with a dataset corresponding to an 
integrated luminosity of ~ 5 fb" 1 of pp collisions collected at y/s = 8 TeV. The LHC8 data 
have shown no evidence for new physics beyond the SM. 

However, there is another side to the story. The concordance model of cosmology -the flat 
expanding universe containing 5% baryons, 20% dark matter, and 75% dark energy- con- 
tinues to be put on a firmer footing through observations of the Supernova Search Team [5- 
7], the Supernova Cosmology Project [8-10], the Wilkinson Microwave Anisotropy Probe 
(WMAP) [11, 12], the Sloan Digital Sky Survey (SDSS) [13-16], and the Hubble Space 
Telescope [17]. While not yet rock solid experimentally, from these observations it is evident 
that in order to describe the physics of the early universe, and thereupon particle interac- 
tions at sub-fermi distances, new theoretical concepts are necessary, which go beyond the 
SM. 1 

Arguably, another major driving force behind the consideration of physics beyond the 
SM is the huge disparity between the strength of gravity and of the SM forces. This hier- 
archy problem suggests that new physics could be at play at the TeV-scale. To be more 
specific, the non-zero vacuum expectation value of the scalar Higgs doublet sets the scale 
of electroweak interactions. However, due to the quadratic sensitivity of the Higgs mass to 
quantum corrections from an aribitrarily high mass scale, with no new physics between the 
energy scale of electroweak unification, Mew ~ 1 TeV, and the vicinity of the Planck mass, 
Mpi ~ 10 19 GeV, the Higgs mass must be fine-tuned to an accuracy of O(10 32 ). Therefore, 
it is of interest to identify univocal footprints that can plausible arise in theories with the 
capacity to describe physics over this colossal range of scales. Among various attempts in 
this direction, string theory is perhaps the most successful candidate and also the most am- 
bitious approach since besides the SM gauge interactions it includes also the gravitational 
force at the quantum level [19, 20]. 

In recent years there has been achieved substantial progress in connecting string theory 
with particle physics and cosmology. Important advances were fueled by the realization of 
the vital role played by D-branes [21, 22] in connecting string theory to phenomenology. 
This has permitted the formulation of string theories with string scale setting in at TeV 
scales, and together with large extra dimensions [23]. 

There are two peerless phenomenological consequences for TeV scale D-brane string com- 
patifications: the emergence of Regge recurrences at parton collision energies \/§ ~ string 
scale = M s ; and the presence of one or more additional U(l) gauge symmetries, beyond 
the U(l)y of the SM. The latter follows from the property that the gauge group for open 
strings terminating on a stack of N identical D-branes is U(N) rather than SU(N) for 
N > 2. (For N = 2 the gauge group can be Sp(l) = SU(2) rather than U(2).) In a series of 
recent publications we have exploited both these properties to explore and anticipate new- 
physics signals that could potentially be revealed at LHC. Regge recurrences most distinctly 
manifest in the 7+ jet [24, 25] and dijet [26-30] spectra resulting from their decay. 2 The 



1 It appears likewise from experimental evidence of neutrino flavor oscillations by the mixing of different 
mass eigenstates that the SM has to be extended [18]. 

2 The amplitudes of lowest massive Regge excitations that include 2 — » 2 scattering processes involving 4 
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extra U(l) gauge symmetries beyond hypercharge have (in general) triangle anomalies, but 
are cancelled by the Green-Schwarz mechanism [33]. In addition there can be also massive 
U(l) gauge bosons, which are associated to 4D non-anomalous Abelian gauge symmetries, 
but however originate from anomalous £/(l)'s in six dimensions. In both cases, these U(l) 
gauge bosons get Stiickelberg masses. Since in these D-brane models M s is assumed to be 
C(TeV), the presence of these generic £7(1) 's may be amenable to experimental tests at the 
LHC [34-36]. 

In this work we take a related but different approach studying new physics effects of 
D-brane models with the conventional assumption TeV <C M s < Mp\. The gauge symmetry 
also arises from a product of U(N) groups, guaranteeing extra U(l) gauge bosons in the 
spectrum. The weak hypercharge is identified with a linear combination of anomalous £/(l)'s 
which itself is anomaly free. As indicated in the preceding paragraph, the extra anomalous 
U(l) gauge bosons generically obtain a string scale Stiickelberg mass. The U(l) symmetries 
survive as global selection rules in the effective low energy theory. Such anomalous gauge 
bosons are now very heavy and out of the LHC reach. 

However, in some D-brane models there exists non-anomalous and also massless U(l) 
gauge symmetries in addition to hypercharge. Namely, under certain topological conditions 
the associated gauge bosons can remain massless and obtain a low mass scale via the ordinary 
Higgs mechanism. Some phenomenological aspects of these kind of U(l) gauge bosons were 
recently discussed in [37]. In this paper we first revisit the prospects of detecting such TeV- 
scale gauge bosons in particular at the LHC, and then we show in all generality that their 
milli-weak interactions could play an important role in observational cosmology. 

Before proceeding with an outline of the paper, we sketch some issues surrounding the 
choice of a non-supersymmetric formulation. To avoid the fine tuning inherent in the hier- 
archy problem, the overwhelmingly favored approach is the introduction of supersymmetry 
(SUSY). However, for the present study, this presents a difficult technical problem: the full 
complexity of the scale of SUSY breaking has been pushed by experiment into the TeV 
region, which coincides with the energy scale involved in searching for the extra U(l) gauge 
bosons. In the absence of an experimental signal for the onset of SUSY breaking, we will 
extract from string theory the choice of the U(l) gauge assignments, as well as the quiver 
structure of the fermionic couplings. In principle, SM-like non-SUSY vacua exist in the 
string landscape [38-40] . Throughout most of this work we will operate within that vacuum 
structure. However, before concluding we will also discuss in some detail the phenomenology 
of supersymmetric vacua and the technical problems associated with a phenomenologically 
viable breaking of an additional U(l) symmetry in a SUSY background. 

The layout of the paper is as follows. In Sec. II we outline the basic setting of intersecting 
D-brane models and discuss general aspects of the effective low energy theory inhereted from 
properties of the overarching string theory. After that, we particularize the discussion to 
the U(3)b x ^(1)^ x U(1)l x U(1)j r intersecting D-brane configuration that realizes the 
SM by open strings [41]. In Sec. Ill we study the associated phenomenological aspects of 



gauge bosons, or 2 gauge bosons and 2 quarks, are universal [26]. Therefore, the s-channcl pole terms of 
the average square amplitudes contributing to 7+ jet and dijet topologies can be obtained independent of 
the details of the compactification scheme. For phenomenological purposes, the poles need to be softened 
to a Brcit-Wigncr form by obtaining and utilizing the correct total widths of the resonances [31]. The 
recent search for such narrow resonances in data collected during the LHC8 run, excludes a string scale 
below 4.69 TeV [32]. 
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non- anomalous U(l) gauge bosons related to experimental searches for new physics at the 
LHC In Sec. IV we explore cosmological predictions of intersecting D-brane models in light 
of recent data, which seem to favor the existence of roughly one additional neutrino species 
(in addition to the 3 contained in the SM), challenging the earliest observationally verified 
landmarks: big bang nucleosynthesis (BBN) and the cosmic microwave background (CMB). 
The gist of Sec. IV extends the previous study of TeV-scale string compactifications [42] 
to D-brane models where some of the U(l) masses are at a high string scale. In Sec. V 
we examine the consequences of possible supersymmetric extensions. Our conclusions are 
collected in Sec. VI. 

II. STANDARD MODEL FROM INTERSECTING D-BRANES 

D-brane string compactifications provide a collection of building block rules that can be 
used to build up the SM or something very close to it [43-57] . In this section, we will briefly 
review the basics of constructing such D-brane models. More comprehensive treatments can 
be found in [58-61]. 

A. Construction Rules and Generalities of D-brane Models 

The details of the D-brane construct depend a lot on whether we use oriented string or 
unoriented string models. The basic unit of gauge invariance for oriented string models is 
a U(l) field, so that a stack of N identical D-branes eventually generates a U(N) theory 
with the associated U (N) gauge group. In the presence of many D-brane types, the gauge 
group becomes a product form JT U(Np), where Np reflects the number of D-branes in each 
stack. As an illustration, consider Type IIA string theory compactified on a six dimensional 
manifold A specific configuration will be given by K stacks of intersecting D6-branes 
filling 4-dimensional Minkowski spacetime M4 and wrapping internal homology 3-cycles of 
*M . Each stack consists of Np coincident D6 branes whose world- volume is M4 x lip, where 
lip is the corresponding homology class of each 3-cycle, with P = 1,...,K. The closed 
string degrees of freedom reside in the entire ten dimensional space, which in addition to 
the gravitational fields, contain the geometric scalar moduli fields of the internal space. 
The open string degrees of freedom give rise to the gauge theory on the D6-brane world- 
volumes, with gauge group Y[U(Np). In addition, there are open string modes which 
split into states with both ends on the same stack of branes as well as those connecting 
different stacks of branes. The latter are particularly interesting: there is a chiral fermion 
living at each four-dimensional intersection of two branes P and Q, transforming in the 
bifundamental representation of U(Np) x U(Nq) [62]. The intersection number of these two 
branes, Ipq = [Hp] ■ [IIq], is a topologically invariant integer whose modulus gives us the 
multiplicity of such massless fermionic content and its sign depends on the chirality of such 
fermions. A particularly simple subfamily of the configurations described above consist of 
taking ^# as a factorizable six-torus: T 6 = T 2 x T 2 x T 2 . We can then further simplify the 
configurations assuming that the 3-cycles can be factorized as three 1-cycles, each of them 
wrapping on a different T 2 . In this case the homology 3-cycle lip can be expressed as 

[Up] = [(np, mp), (np, m 2 P ), (np, mp)] , (1) 
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where (rip, m P ), are the wrapping numbers of each D6p-brane, on the i th torus, with n P and 
(m P ) being the number of times the brane is wrapping around the i th torus. The intersection 
number takes a simple form 

3 

!pq = II( n ^ m Q _ ml p n Q) ■ ( 2 ) 
i=i 

In orientifold brane configurations, which are necessary for tadpole cancellation [63, 64], 
and thus consistency of the theory, open strings become in general non-oriented. For unori- 
ented strings the above rules still apply, but we are allowed many more choices because the 
branes come in two different types. There are the branes whose images under the orientifold 
are different from themselves and their image branes, and also branes who are their own 
images under the orientifold procedure. Stacks of the first type combine with their mirrors 
and give rise to U (N) gauge groups, while stacks of the second type give rise to only SO(N) 
or Sp(N) gauge groups. 

Generally speaking, intersecting D-brane models involve at least three kinds of generic 
mass scales. First, of course, there is the fundamental string scale, 

M. = -L, (3) 



where a' is the slope parameter of the well known Regge trajectories of vibrating strings 

j=j + a'M\ (4) 

with j and M = y/nM s the spin and mass of the resonant state, respectively (n = 1, . . . ). 
Second, compactification from ten to four dimensions on an internal six-dimensional space 
of volume Vq defines a mass scale: 

M 6 = -^- (5) 

Third, wrapping a stack P of D(p + 3)-branes around the internal p-cycle defines an internal 
world-volume Vp P ^ = (27r) p Vp of this D-branes stack and an associated (Kaluza-Klein) 
mass: 

MP = i-n- . (6) 

(4 P| ) 

These three types of fundamental dimensional parameters of D-brane models are linked 
to four-dimensional physical observables. First, the Planck mass given by 

Ml, = 8 Ml ^ (7) 

determines the strength of gravitational interactions. Here, the dilaton field 0io is related 
to the string coupling constant through g s = e^ 10 . Thus, for a string scale M s m (9(1 TeV), 
the volume of the internal space M 6 needs to be as large as V^M® = C(10 32 ). Second, the 
four-dimensional gauge couplings of the strong and weak interactions are given in terms of 
the respective volumes V p p , where P runs over the corresponding gauge group factors, as 

g P 2 = (27T)" 1 MP . (8) 
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Again for a string scale M s TeV) and using the known values of the strong (g 2 /47r 

0.1) and the weak (<? 2 /47r « gf/127r) gauge coupling constants at the string scale (grf / 47r = 
qem/ sin 2 9w, sin 2 9w ~ 0.23, «em ~ 1/128) we can compute the volumes of the internal 
cycles, assuming weak string coupling. To be specific, we choose g s = 0.2, and then we 
obtain 

Mf « 1 , M? v® « 3 . (9) 

For TeV <C M s < Mpi, V(j and Vp^'s are 0(1) in string units. In general, there are different 
volumes V p s for different stacks, and therefore the abelian gauge couplings associated to 
£7(1) symmetries of different D-brane stacks are not equal. 

This approach to string model building leads to a variety of low energy theories including 
the SM as well as its supersymmetric extensions. Throughout most of this paper we consider 
theories which are non-supersymmetric all the way up to the UV cutoff of the effective theory; 
of course the deep UV theory of quantum gravity may well be supersymmetric. Even though 
SUSY introduces special advantages over completely non-SUSY theories, our approach is 
distiguished by its simplicity to describe very appealing phenomenological possibilities that 
best display the dynamics involving the extra U(l) symmetries. The study of some aspects 
of the supersymmetric version of these models will be postponed until Sec. V. 

The minimal embedding of the SM particle spectrum requires at least three brane 
stacks [65] leading to three distinct models of the type £7(3) x £7(2) x £7(1) that were clas- 
sified in [65, 66]. Only one of them (model C of [66]) has baryon number as symmetry 
that guarantees proton stability (in perturbation theory), and can be used in the framework 
of TeV-scale strings. Moreover, since the charge associated to the £7(1) of U(2) does not 
participate in the hypercharge combination, 17(2) can be replaced by the symplectic Sp(l) 
representation of Weinberg-Salam SU(2)l, leading to a model with one extra £7(1) (the 
baryon number) besides hypercharge [67]. 

The SM embedding in four D-brane stacks leads to many more models that have been 
classified in [68, 69]. In order to make a phenomenologically interesting choice, herein we 
focus on models where £7(2) can be reduced to Sp(l). Besides the fact that this reduces the 
number of extra £7(1) 's, one avoids the presence of a problematic Peccei-Quinn symmetry, 
associated in general with the £7(1) of U(2) under which Higgs doublets are charged [65]. 
To develop our program in the simplest way, we will work within the construct of a minimal 
model, U(3)b X Sp(1)l X U(\)l x U(1)j r , which has the attractive property of elevating the 
two major global symmetries of the SM (baryon number B and lepton number L) to local 
gauge symmetries [41]. We turn now to discuss the compelling properties of this model. 

B. Standard Model++ 

In this paper we are interested in the minimal 4-stack gauge-extended sector U(3)b X 
Sp(1)l xU(l)iX U(l)i R [41]. A schematic representation of the D-brane structure is shown 
in Fig. 1 and the brane content is given in Table I. Note that for the Sp(l) stack P, the 
mirror brane P* lies on top of P. So even though Np = 1, it can be thought of as a stack 
of two D6 branes, which give an Sp(l) = SU(2) group under the orientifold projection. 
Concretely, in the bosonic sector the open strings terminating on the stack of "color" branes 
contain, in addition to the SU(3) octet of gluons 

G% = [pfjGl - d u Gl + g 3 f abc G b ^G c u ) , if abc T a = [T\ T c ] , T a e SU(3) , 
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Sp(D U(l) 

FIG. 1: Pictorial representation of the U(3)b x Sp(1)l x U(1)l x U(1)i r D-brane model. 

an extra U(l) boson C M . On the S'p(l) stack the open strings correspond to the weak gauge 
bosons 

W; u = faW? - d v W£ + ^e^WjW?) , ie abc r a = [t\ t c ] , r a = a a /2 e SU(2) . 

The U(l)i a D-brane is a terminus for the B^ gauge boson, and there is a third additional 
U(l) field terminating on the U(1)l brane. The resulting U(l) content gauges B [with 
U{1)b C U(3)b], L, and a third additional abelian charge Ir which acts as the third isospin 
component of an SU(2)r. The usual electroweak hypercharge is a linear combination of 
these three U(l) charges: 

Qy = c\Qi R + c 3 Q B + CiQ L , (10) 

with ci = 1/2, c 3 = 1/6, C4 = —1/2, B = Qb/3 and L = Ql. Alternatively, inverting the 
above relations, one finds: 

Q B = 3B ; Q L = L ; Q Ir = 2Q Y - (B - L) . (11) 

The chiral particle spectrum from these intersecting branes consists of six sets (labeled by 
an index % — 1, . . . , 6) of Weyl fermion-antifermion pairs, whose quantum numbers are given 
in Table II. Note that the combination B — L is anomaly free, while both B and L are 
anomalous. 

As mentioned already, the Qb (gauged baryon number) is anomalous. This anomaly is 
canceled by the 4D version [70-74] of the Green-Schwarz mechanism [33]. Non anomalous 
£7(1) 's can acquire masses due to effective six- dimensional anomalies associated for instance 
to sectors preserving M = 2 supersymmetry [75, 76]. 3 These two-dimensional 'bulk' masses 



In fact, also the hypercharge gauge boson of U(l )y can acquire a mass through this mechanism. In order 
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TABLE I: D-brane content of U(3) B x Sp(l) L x U(1) L x C/(l)/ H ; the mirror branes 0*,P*,Q*,R* 
are not shown. 



Label 


Stack 


Number of Branes 


Gauge Group 


1 = R 


Right 


N R = 1 


U(l)i R 


2 = P 


Left 


N P = 1 


Sp(l) L * SU(2) L 


3 = Q 


Baryonic 


Nq = 3 


U(3) B = SU(3) C x U(1) B 


4 = 


Leptonic 


N = l 


U(1)l 



TABLE II: Chiral fermion spectrum of the U(3) B x Sp(1)l X U(1)l x U(1)i r D-brane model. 



Label 


Fields 


Sector 


IPQ 


Representation 


Qb 


Ql 


QlR 


Qy 


1 


Ur 


(Q,R*) 


3 


(3,1) 


1 





1 


2 
3 


2 


Dr 


(Q,R) 


3 


(3,1) 


1 





-1 


1 

3 


3 


Ll 


(o,p) 


3 


(1,2) 





1 





1 

2 


4 


Er 


(0,R) 


3 


(1,1) 





1 


-1 


-1 


5 


Ql 


(Q,P) 


3 


(3,2) 


1 








1 
6 


6 


Nr 


(0,R*) 


3 


(1,1) 





1 


1 






become therefore larger than the localized masses associated to four- dimensional anomalies, 
in the large volume limit of the two extra dimensions. Specifically for D(p + 3)-branes 
with p-longitudinal compact dimensions the masses of the anomalous and, respectively, the 
non-anomalous U(l) gauge bosons have the following generic scale behavior: 

anomalous U(l) a : M z > = g' a M s , (12) 
non — anomalous U(l) a : M Z " = g' a Mg V 2 . 

Here g' a is the gauge coupling constant associated to the group U (l) a , given by g' a oc g s / \Jvj^ 

where g s is the string coupling and Vp a ^ is the internal D-brane world-volume along the p 
compact extra dimensions, up to an order one proportionality constant. Moreover, V% is 
the internal two-dimensional volume associated to the effective six- dimensional anomalies 
giving mass to the non-anomalous t/(l) a . 4 E.g. for the case of D5-branes, whose common 
intersection locus is just 4-dimensional Minkowski-space, Vp = V% denotes the volume of 



to keep it massless, certain topological constraints on the compact space have to be met. 
4 It should be noted that in spite of the proportionality of the U(\) a masses to the string scale, these 
are not string excitations but zero modes. The proportionality to the string scale appears because the 
mass is generated from anomalies, via an analog of the Green-Schwarz anomaly cancellations: cither 4 
dimensional anomalies, in which case the Green-Schwarz term is equivalent to a Stiickelberg mechanism, 



s 



the longitudinal, two-dimensional space along the two internal D5-brane directions. Since 
internal volumes are bigger than one in string units to have effective field theory description, 
the masses of non-anomalous ?7(l)-gauge bosons are generically larger than the masses of 
the anomalous gauge bosons. 5 

The non-anomalous U(l) a can also remain massless all the way down to the TeV-scale 
energy region and grow a mass through a Higgs mechanism. The absence of a Stiickelberg 
mass term for the associated gauge bosons means that this U(l) gauge symmetry is anomaly 
free also in six dimensions. In this case a certain topological condition has to hold, which 
cannot be read off from the local D-brane quiver, but can only be answered knowing the 6D 
compact orientifold. Specifically, just like for the SM gauge symmetry U(1)y, the absence 
of the Stiickelberg mass term for U(l) a = c%U(1)b + c%U{1)l + c"t/(l)/ fl can be phrased 
by the following condition on the homology cycles II and their orientifold images II' of the 
three U(l) gauge groups: 

3c^(n 3 - n' 3 ) + c a 4 (u L - u' L ) + ct(u lR - ny = o . (is) 

In what follows we entertain this possibility, having two massless gauge bosons U(1)y 
(associated to the SM hypercharge) and U(1)y" (associated to a linear combination of 
anomaly-free Ir and B — L) and one heavy gauge boson U(1)y' (associated to an anomalous 
combination of the three ?7(l)'s). The classical gauge invariant Lagrangian, obeying the 
U(3)b x Sp(1)l x U{1)l x U{l)i R gauge symmetry, can be decomposed as: 

•^sm ++ = =^ym + (~^f °^ Y ) + + , (14) 

generations 

where the terms on the right hand side identify the gauge (or Yang-Mills) part, the fermion 
part, the Yukawa part, the scalar part, and extra terms from the underlying string theory, 
respectively. 

Electroweak symmetry breaking is achieved through the standard Higgs doublet H. The 
spontaneous symmetry breaking of the extra non-anomalous U(l) is attained through an 
SU (2) singlet scalar field H", which acquires a vacuum expectation value (VEV) at the TeV 
scale. The U(l) quantum numbers of the Higgs sector are given in Table III. 

The Yang-Mills Lagrangian reads: 

^ym = -\ (G^GT + W^WjT + F$F$ + F$Ffi + FWFgj) , (15) 

with the non-Abelian field strengths the same as in the SM, and the Abelian F^J = d^B v - 

d v B^ F$ = dpC v - d v C„ and F$ = d^X u - d u X^. 
The fermion Lagrangian is given by 

&f = iQh^Qi + iUhvP^Un + iD~ Rl p^D R + i~T Ll ^ L L + iE^y^En 

+ iNr^Nr, (16) 



or from effective 6 dimensional anomalies, in which case the mass term is extended in two more (internal) 
dimensions. 

5 In [77] a different (possibly T-dual) scenario with D7-branes was investigated. In this case the masses of 
the anomalous and non-anomalous U(l)'s appear to exhibit a dependence on the entire six-dimensional 
volume, such that the non-anomalous masses become lighter than the anomalous ones. 
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TABLE III: Higgs spectrum of the U{3) B x Sp(l) L x U{l) L x U(1) Ir D-brane model. 



Fields 


Sector 


IPQ 


Representation 


Qb 


Ql 


QlR 


Qy 


H 


(P,R) 


1 


(1,2) 








1 


1 

2 


H" 


(0,R) 


1 


(1,1) 





-1 


-1 







where 

V, = d, - ig 3 T a G; - ig' 3 Q B C, - tg 2 r a W; - ig' x Q lR B^ - ig' 4 Q L X, (17) 

are the covariant derivatives with the gauge fields specified in the D-brane basis. 
The fields C^, X^, are related to Y^, and YJ' by the rotation matrix, 

C^C^ + S^SgS^p —S^C^ + C^SgS^ I , (18) 
S^Cg C^Cg } 

with Euler angles 9, ip, and [78]. Hence, the covariant derivative for the U(l) fields in 
Eq. (17) can be rewritten in terms of Y" M , Y', and Y" as follows 

V, = 8,- z% (-S e g[Qi R + C e S^g' 4 Q L + CgC,pg' 3 Q B ) 

~ [CoScj>9iQi R + (CpC^, + SgSpS^,) g' 4 QL + (C^SgSp — C^S^g^Qs] (19) 
— iY^ [CgC ( f > g l Qi R + (—CipS<p + C^SgS^p) g' 4 Qh + (C^C^Sg + S^S^) g' 3 Q B ] ■ 

Now, by demanding that Y^ has the hypercharge Qy given in Eq. (10) we fix the first column 
of the rotation matrix IR 

Y^c 3 g Y /g' 3 • • • \ 

Y,c 4 g Y /g' 4 ... , (20) 
y^cig Y /g[ ■■■ J 
and we determine the value of the two associated Euler angles 

9 = -arcsinfcigy/gi] (21) 

and 

if) = arcsin[c 4 sfy/ (g' 4 C e )} . (22) 

The couplings g[ and g' 4 are related through the orthogonality condition, P(g Y , g[, g'3, g'4) = 0, 
yielding 

s)H-(t) 2 -(t) 2 ' 

with g' 3 fixed by the relation g^(M s ) = y/Qg' 3 (M s ) [36]. Next, by demanding that Y" couples 
to a linear combination of anomaly-free I B and B — L we determine the third Euler angle 

tan0 = -S e v v . (24) 

6 g 3 b^p — g 4 U^f, 

10 



In the (Y, Y', Y") basis, Y and Y" are coupled to anomaly-free currents while the anomaly 
of the current associated to Y' is cancelled by the generalized Green-Schwarz mechanism. 
As a result, Y' acquires a mass of order of the string mass M s , c.f. Eq. (12). Higgs VEVs 
will generate additional mass terms for Y', introducing also some small mixing with other 
gauge gauge bosons, of order (TeV/M s ) 2 . From now on, we neglect such small effects and 
take Y' ~ Z' . 

The Yukawa interactions are given by 

= ~Y d (QlH) D r - Y u (Q2ta 2 H*) U R - Y e (T L H) E R - Y N (LT L ia 2 H*) N R + h.c, (25) 

where the Yukawa couplings Yi are matrices in flavor space. Note that unlike in the super- 
symmetric case, a single Higgs vacuum expectation value will generate masses for up and 
down quarks. 6 

Note that with the charge assignments of Tables II and III there are no dimension 4 
operators involving H" that contribute to the Yukawa Lagrangian. This is very important 
since H" carries the quantum numbers of right-handed neutrino and its VEV breaks lepton 
number. However, this breaking can affect only higher-dimensional operators which are 
suppressed by the high string scale, and thus there is no phenomenological problem with 
experimental constraints for M s higher than ~ 10 14 GeV. 

The scalar Lagrangian is 



= (V^HY V^H + {V^H") ] V^H" - V(H, H") , 

with the potential 

V (H, H") = /i 2 \H\ 2 + // 2 \H"\ 2 + Ai \H\" + A 2 \H"\ 4 + A 3 \H\ 2 \H"\ 2 
The Higgs VEVs obtained after minimizing this potential will be denoted as 



(26) 



(27) 



(H) 



and (H") 



(28) 



The kinetic terms of the Higgs fields in (26) give masses to the various gauge bosons. 

At this point, we identify the photon and weak force mediators W+, W~, perform- 
ing the usual Weinberg rotation 



( Cq w Sg w 

-Sg w Cg w 



\ (Y,\ 



1/V2 i/V2 

1/V2 -1/V2) 



this gives 



W 3 



(29) 



V 



X 



d hl - - 92 °~W+ - -92 v + W~ - 192 cos 9 W (a 3 /2 - Q Y tan 2 9 W ) Z ^ - ig 2 sin W 



(a 3 / 2 + Qy) Af, - igY'Qv'Z'^ - ig Y »QY»Y{{ , 



(30) 



i<j2H* transforms in the fundamental representation of SU(2). 
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with = (cr 1 ± ia 2 ) /2 , gy/g2 = t&nOw- From (19) and (30) we define 

Q Y H = if/2, 

gyiQyiH = (g^gS^H , 

gynQynH = {g^CeC^H , 

Qy#" = 0, 

9y'Qy'H" = —[g'lCgStp + g'^C^Cif, + SgS^S^)]!!" , 

gynQynH" = —(g[CeC^, + g'^C^SgSjj, — C^Sfj,)]!!" . 



(31) 



The Higgs kinetic terms of Eq.(26) together with the Green-Schwarz mass term, \M' 2 Z'^ , 
lead to 



= [Z>], (0 «)] 



(?V/')t(W) 



M' 2 Z[ Z'v 



(32) 



Expanded this gives 
1 



+ v" 2 {g[C e (S^ Z; + + ^ [(C^ + SgS^) % + S^C^ F^] } 



* 4G& v) 2 W+W~» + -^vfC^ Z„Z* + g[CgC^g 2 v 2 C^Z 
+ v" 2 (g[C e C^ Y; + g^S e C^ Y^f + (g[v C CtfY'p"* + ... 



(33) 



where the omitted terms pertain only to the Z' couplings at the string scale. Recall that we 
have taken M' ~ M s and therefore Z' decouples from the low energy physics. By inspection 
of (33) we immediately recognize the W ± masses and the usual tree level formula for the 
mass of the Z particle in the electroweak theory, M z = (g 2 v 2 + g Y v 2 )/2, before mixing. 
Now, we use the relation g[Sg = g'^CgS^ to conveniently rewrite (33) as 



a ^ \ (92V) 2 w+w-* + \{ 9 2v) 2 c~ 2 zX + g [vc e c, 92 vc^ y>;z» + ( 



y"g[CA 
c. 



{g 2 v) 2 W+W-» 



G 



+ 



-{g 2 vfC7 2 - 9[ Cl °\ 92 ^ °^ 

4 w HV'gnf \i + &Cj) 




g[ CI C+ g 2 v 2 C£ Z, 



2(v"g' 1 C 4> ) 2 



Z u Z» + ... . 



(34) 



Finally, if we make the expansion around v/v" <C 1, the Z^Y'" 1 mass matrix is render 
diagonal and we obtain the desired expression for the mass terms 



g2v 



2C, 



(35) 
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where Z" ~ Y" + small corrections. 

In principle, in addition to the orthogonal field mixing induced by identifying anoma- 
lous and non-anomalous U(l) sectors, there may be kinetic mixing between these sectors. 
However, in models where there is only one U{1) per stack of D-branes, the relevant kinetic 
mixing is between C/(l)'s on different stacks, and hence involves loops with fermions at brane 
intersection. Such loop terms are typically down by gf/16ir 2 ~ 0.01 [79]. 7 By inspection 
of Table II the charges Qb, Ql, and Qj R are mutually orthogonal in the fermion space, 
i.e QijQjj = for i ^ j. This will maintain the othogonality relation P = to one 
loop without inducing kinetic mixing [36]. The charges assigned to H" (see Table III) will 
violate the orthogonality condition. However, the H" only contributes at the 0.9% level to 
the running of g[ from the string scale to the TeV scale, and about 0.3% to the running of 
g 4 . These are of the same order as the two loop contributions from the fermion sector, so we 
may consistently ignore the nonorthogonality introduced by H" in the context of one loop 
considerations. 

III. LHC PHENOMENOLOGY 

In this section we discuss the discovery potential of the Z" resonance at the LHC. Before 
proceeding, we summarize the lessons learned thus far. The initially free parameters of the 
model consist of three couplings g[, g' 3 , g' 4 . These are augmented by three Euler angles to 
allow for a field rotation to coupling diagonal in hypercharge. This diagonalization fixes two 
of the angles and the orthogonal nature of the rotation introduces one constraint on the 
couplings P(gy, 9i, g'3, gi) — 0. The baryon number coupling g' 3 is fixed to be ^/l/G of the 
non-abelian 577(3) coupling at the scale of U(N) unification, and is therefore determined at 
all energies through RG running. In what follows, we take M s = 10 14 GeV as a reference 
point for running down the g' 3 coupling to the TeV region that is ignoring mass threshold 
effects of stringy states. This yields g' 3 (M s ) = 0.231. We have checked that the running of 
the g' 3 coupling does not change significantly within the LHC range, for different values of the 
string scale. This leaves one free angle and two couplings with one constraint. Equation (24) 
fixes the third Euler angle. To comply with these assignments and ensure perturbativity of 
g' 4 between the TeV scale and the string scale we find from (23) that g[ > 0.4845. We also 
take g[ < 1 in order to ensure perturbativity at the string scale. 

We first consider the case with g[(M s ) ~ 1. This leads to ip(M s ) = -1.245, 6(M S ) = 
-0.217, <f){M s ) = -0.0006, and g 4 (M s ) = 0.232. Substituting our fiducial values in (19) we 
find the non-anomalous U(l) vector bosons couple to currents 

J y = 2.1 x 10" 1 Q Ir + 2.1 x 10" 1 (B - L) 

J Y „ = 9.8 x 10" 1 Q Ir - 4.7 x 10~ 2 (B-L), (36) 

at the string scale. Next, we run the couplings down to the TeV region. A very important 
point is that the couplings that are running are those of the U(l) fields; hence the (3 functions 
receive contributions from fermions and scalars, but not from gauge bosons. The one loop 



7 The major effect of the kinetic mixing is in communicating SUSY breaking from a hidden U(l) sector to the 
visible sector, generally in modification of soft scalar masses. For a comprehensive review of experimental 
limits on the mixing, see [80]. 
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correction to the various couplings are 



^ HQ /Ms) , (37) 



a Y (Q) a Y (M s ) 2tt 

- HQ/Ms) , (38) 



ai(Q) a,(M s ) 2tt 
where 

with / and s indicating contribution from fermion and scalar loops, respectively. Setting 
Q = 4 TeV, from (38) we obtain: g[ = 0.406, g' 3 = 0.196, ^ = 0.218, 9 = -0.466, 
V> = -1.215, and <p = -0.0003. This leads to 

J Y = 1.8 x 10" 1 Q Ir + 1.8 x 10" 1 (B - L) 

J Z " = 3.6 x lO' 1 Q Ir - 9.2 x 10" 2 (B - L) , (40) 

where we have assumed that H" has developed its VEV. Since Tr [Qi R B] = Tr [Qj R L] = 0, 
the Z" decay width is given by 

oc (1.4 x 10" 1 ) 2 Tr[Qj R } + (9.2 x 10~ 2 ) 2 Tr [(B - L) 2 ] 

= 1.0 x 10° + 4.5 x 10" 2 . (41) 

Thus, the corresponding branching fractions are BR Z" — > Qj R = 0.959 and BR Z" — > 
B — L = 0.041. Though not relevant for LHC phenomenology, a straightforward calculation 
shows that Z' is very nearly diagonal in B, with BR Z' — >■ B = 0.946 and BR Z' — > L = 
0.054. Of course, since the quiver construction has each particle straddling two adjacent 
branes, there can be considerable variation in decay channels particle by particle. This is 
evident in Table IV. 8 The dominance of B for the Z 1 decay channel and Ir for the Z" decay 
channel is valid after averaging over decay channels. 9 

Now, duplicating the procedure for g[(M s ) = 0.4845 we obtain 

BR Z' -> B : BR Z' -> L : BR Z" -»■ Q m : BR Z" -+ B - L . , 

0.066 : 0.934 : 0.039 : 0.961 . ^ ' 

The chiral couplings of Z' and Z" gauge bosons which are mostly L and B — L, respectively 
are given in Table V. The variation of the branching fractions within the allowed range of 
g[(M s ) is shown in Fig. 2. 



The physical couplings of the Z" to fermions helds given in Table IV are consistent with the bounds 
presented in [81] from a variety of experimental constraints. 

An analogue is in the SM. The Z couples to a current Jz oc I3 — tan 2 dw\, where Q = T3 — In this 
case, E(t) 2 = T and Tr p3 ] = 2 5 we have BR z ~> T 3 ■ BR Z -> | = 2 : f tan 4 6> w = 2 : 0.25 = 8:1. 
However, this certainly does not hold particle by particle; e.g., for the neutrino electron doublet: r z^v oc 
(1 + tan 2 6* w ) 2 ~ 1.7, whereas T z ^ e oc (1 - tan 2 6 W ) 2 - 0.5. 
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FIG. 2: Branching fractions of Z' (left) and Z" 
denote the branching into B (left) and Ir (right) 
(left) and B-L (right). 
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(right) as a function of g[(M s ). The solid lines 
. The dashed lines denote the branching into L 



The LHC discovery potential for Z" gauge boson as a mass peak above a small background 
in the reactions pp — > Z" — >■ jj and pp — > Z" — > with I = e,p, is well known. The 

required luminosity to discover a Z" basically depends only on its cross section, and therefore 
on its mass and couplings. Experimental effects due to mass resolution are known to result 
in an only minor reduction of the sensitivity. 

Using a data set of pp collisions at \fs = 8 TeV, with an integrated luminosity of 4.0 fb _1 , 
the CMS Collaboration has searched for narrow resonances in the dijet invariant mass spec- 
trum [32]. Each event in the search is required to have its two highest-pT jets with (pseu- 
dorapidity) \r)j\ < 2.5. The acceptance A of selection requirements is reported to be ~ 0.6. 
The spectra are consistent with SM expectations and thus upper limits on the cross section 
times branching fraction for Z" into two jets have been set. Similar upper limits have been 
obtained by the ATLAS Collaboration using 5.8 ib' 1 of data collected at yfi = 8 TeV [82]. 
These results, which are display in Fig. 3, extend previous exclusion limits from LHC7 [83- 
87]. 

The ATLAS Collaboration has searched for narrow resonances in the invariant mass 
spectrum of dimuon and dielectron final states in event samples at = 7 TeV corresponding 
to an integrated luminosity of 4.9 fb" 1 and 5.0 fb" 1 , respectively [88]. The spectra are 
consistent with SM expectations and thus upper limits on the cross section times branching 
fraction for Z" into lepton pairs have been set. More recently, the CMS Collaboration 
updated the LHC7 results using 4.1 fb _1 of data collected at y/s = 8 TeV [89]. The combined 
upper limits from LHC7 and LHC8 are shown in Fig. 4. Previous dilepton searches by the 
LHC experiments have been reported in [90, 91]. 

In order to compare with these results we now turn to compute these cross sections in 
our model. The relevant part of (16), ffZ" coupling, is of the form 

= E flYfl + (9y»Qy») fk fkffk) Z'> , (43) 

where /^(R) are f erm i° n fields and £f Lt f R = v q ± a q , with v q and a q , the vector and axial 
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TABLE IV: Chiral couplings of Y, Z', and Z" gauge bosons. All fields in a given set have a 
common qy'Qy' ■> 9y"Qy" couplings. We have taken Z' to be mostly B and and Z" to be mostly 
Ir- 



Fields 


9yQy 


qy'Qy 1 


9y"Qy" 


Ur 


0.2434 


0.1836 


0.3321 


Dr 


-0.1214 


0.1838 


-0.3933 


U 


-0.1826 


0.0759 


0.0918 


Er 


-0.3650 


0.0760 


-0.2709 


Ql 


0.0610 


0.1837 


-0.0306 


Nr 


0.0000 


0.0758 


0.4545 


H 


0.1824 


0.0000 


0.3627 


H" 


0.0000 


-0.0758 


-0.4545 



TABLE V: Chiral couplings of Y, Z', and Z" gauge bosons. All fields in a given set have a common 
9Y'Qy'i 9y"Qy" couplings. We have taken Z' to be mostly L and and Z" to be mostly B — L. 



Fields 


9yQy 


9y'Qy> 


9y"Qy" 


Ui 


0.2435 


0.1101 


-0.0763 


Di 


-0.1217 


0.1101 


-0.2242 


Li 


-0.1825 


0.7165 


0.4509 


Ei 


-0.3651 


0.7165 


0.3769 


Qi 


0.0609 


0.1101 


-0.1503 


N, 


0.0000 


0.7165 


0.5248 


H 


0.1826 


-0.0000 


0.0739 


H" 


-0.0000 


-0.7165 


-0.5248 



couplings respectively. To compare our predictions with LHC experimental searches in 
dilepton and dijets it is sufficient to consider the production cross section in the narrow Z" 
width approximation, 

a(qq -+ Z") = I<f [v^gi) + a^g',)] 6 (s - M*,,) , (44) 

where Gf is the Fermi coupling constant and the A'-factor represents the enhancement from 
higher order QCD processes estimated to be K ~ 1.3 [92]. After folding a with the CTEQ6 
parton distribution functions [93], we determine (at the parton level) the resonant production 
cross section. In Figs. 3 and 4 we compare the predicted cr(pp — > Z") x BR(Z" — > jj) and 
a(pp — > Z") x BK(Z" — > ££) production rates with 95% CL upper limits recently reported 
by the CMS and ATLAS collaborations. Selection cuts will probably reduce event rates by 
factors of 20%. Keeping this in mind, we conclude that if Z" is mostly Ir, then the predicted 
production rates for M Z " ~ 4 TeV at yfs = 8 TeV saturate the current dijet limits. On the 
other hand, if Z" is mostly B — L the lower limit on the gauge boson mass, Mz» ^ 3 TeV, 
is determined primarily from dilepton searches. 

For the discovery potential in the high mass region the dijet channel is statistically a 
better discriminator than lepton pairs. Therefore, we investigate (at the parton level) the 



1(3 




FIG. 3: Comparison of the (pre-cut) total cross section for the production of pp — > Z" — > jj with 
the 95% CL upper limits on the production of a gauge boson decaying into two jets as reported by 
the CMS and ATLAS collaborations (corrected by acceptance). For isotropic decays (independently 
of the resonance), the acceptance for the CMS detector has been reporetd to be A ~ 0.6. The 
ATLAS acceptance ranges from 11% to 54% varying from 1 TeV to 4.25 TeV, and is never lower 
than 48% for masses above 2 TeV. The case in which Z" is mostly diagonal in Jr is shown in the 
left panel and the case in which it is mostly B — L in the right panel. 




FIG. 4: Comparison of the (pre-cut) total cross section for the production of pp —> Z" 11 with 
the 95% CL upper limits on the production of a gauge boson decaying into two leptons, as reported 
by the ATLAS and CMS collaborations. The case in which Z" is mostly diagonal in Ir is shown 
in the left panel and the case in which it is mostly B — L in the right panel. 
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LHC14 sensitivity for a Z" resonance (which is mostly Ir) in data binned according to the 
dijet invariant mass M, after setting cuts on the different jet rapidities, \yx\, 1 3/32 ] < 1 and 
transverse momenta p^. 2 > 50 GeV. With the definitions Y = \{yi + 2/2) and y = \{yi — 2/2), 
the cross section per interval of M for pp — > dijet is given by 

- = Mr 

ijkl 



dM 



Mr I dY M x a, M) fj(x b , M) I 

AAU1 — ^rnax ** 



dy 



^raax r ^max ^ 

dY fi(x a ,M) fj(x b ,M) / dy 



(j/max + V) 

da 



da 
dt 

1 



ij^ki cosh y 



d cosh y 



(45) 



where f(x, M)'s are parton distribution functions (we use CTEQ6 [93]), r = M 2 /s, x a 
\fre Y , = y/re~ Y , and 



|A4(z'j -> A;/)| 2 = 16tts 2 



da 
di 



(46) 



we specify partonic subprocesses with caret notation (s, t, u). The y integration range in 
Eq. (45), Fmax = min{ln(l/ v /r), 2/ max }, comes from requiring x a , x b < 1 together with the 
rapidity cuts y min < |yi|, |y 2 | < Umax- The kinematics of the scattering also provides the 
relation M = 2ppcosh?/, which when combined with p? = M/2 sin#* = M/2y/l — cos 2 6*, 
yields coshy = (1 — cos 2 #*) -1 / 2 , where 9* is the center-of-mass scattering angle. Fi- 



nally, the Mandelstam invariants occurring in the cross section are given by s 



M 2 



t 



■\M 2 e 



1 / cosh?/, and u = —\M 2 e +y / cosh?/. 



The average square amplitude (for incoming quark q and outgoing quark q') is given by, 



\M(qq^q'q')\ 2 



X 



- [g Y "Q Y "(QL) + 9y"Qy"(<1r)] [9y"Qy"(Ql) 
2(u 2 + t 2 ) 



g Y »QY»( a R)] 



M 2 „) 2 + (T Z n M z „) 2 



(47) 



where gY"QY"{.<iL) an d 9y"Qy"{.<1r) are the couplings of Z" to quarks. (Note that we have 
not summed over the flavors, but we did average and sum the colors). 
The decay width of Z" f f is given by 



r(Z" -+ //) = <f^N e Mz»VT=te [vj(l + 2x) + aj(l - Ax)] (l + -) , 



6ttV2 



(48) 



where ot s = a s (Mz>>) is the strong coupling constant at the scale Mz», x = m 2 /M^„, Vf 
and a/ are the vector and axial couplings, and N c = 3 or 1 if / is a quark or a lepton, 
respectively [94]. For our fiducial values of g[ and <p we obtain v 2 + a 2 = 0.396 and v 2 d + a 2 , = 
0.554. 

We calculate a signal-to-noise ratio, with the signal rate (S) estimated in the invariant 
mass window [M^" — 2T, Mz» +2F]. To accommodate the minimal acceptance cuts on dijets 
from the CMS and ATLAS proposals [95], an additional kinematic cut, \y max \ < 1-0, has 
been included in the calculation. The noise (AT) is defined as the square root of the number 
of QCD background events (B) in the same dijet mass interval for the same integrated 
luminosity. In Table. VI we show the behavior of the signal-to-noise ratio as a function of 
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the mass of Z" at LHC14 for different integrated luminosities. We conclude that the LHC 
provides a generous discovery potential for Z" which is mostly Ir. The discovery potential 
of a Z" which is mostly B — L is controlled by the sensitivity of LHC14 to dilepton final 
states. For 300 fb _1 , the projected sensitivity is M Z " < 5 TeV [96]. 



TABLE VI: Signal-to-noise ratio at LHC14 for different integrated luminosities. 







10 fb' 


l 




100 ftr 


-l 




iooo ftr 


i 


M z „ (TeV) 


S 


B 


S/N 


S 


B 


S/N 


S 


B 


S/N 


4 


39 


579 


1.62 


391 


5789 


5.14 


3910 


57895 


16.25 


5 


7 


176 


0.50 


67 


1759 


1.60 


670 


17590 


5.05 


6 


1 


66 


0.14 


11 


664 


0.44 


113 


6646 


1.39 



If the Z" is observed at the LHC, we will obviously measure its mass, its total width and 
cross section. In addition, the off- and on resonance peak forward-backward charge asymme- 
tries Ap B would provide additional information about Z" couplings and interference effects 
with the Z boson and the photon [97-99]. 10 Besides, the Z" rapidity distribution is sensitive 
to the uuZ" and bbZ" gauge couplings. Since the W ± and Z boson rapidity distributions 
will be measured in great detail at the LHC, rapidity spectra for the mass region of interest 
can be calculated separately for uu, dd, and sea quark antiquark annihilation. A combined 
fit to the relative parton distribution functions and the Z" rapidity distribution would allow 
us to obtain the fraction of Z" bosons produced from uu and dd initial states [100]. 

IV. NEUTRINO COSMOLOGY REDUX 

In this section we reexamine some critical cosmological issues surrounding the presence 
of the six additional neutrino degrees of freedom correlated to the presence of Z" in our 
dynamical D-brane model. These considerations, when viewed in the context of most re- 
cent cosmological observations are found to constrain the mass of the Z" to an interesting 
narrow band, which will be directly probed by LHC14. To provide a starting point, we first 
summarize the "best-fit" cosmological parameters to recent data. 

A. Beyond ACDM 

Our universe seems, according to the present-day evidence, to be spatially flat and to 
posses a non- vanishing cosmological constant (A) plus cold dark matter (CDM), corre- 
sponding respectively to roughly 70% and 25% of the total density, with the remaining 
5% in baryons. The standard ACDM cosmology provides a rather good fit of existing data 
from BBN (~ 20 minutes), the CMB (~ 380 kyr), and the galaxy formation epochs of the 



The leptonic forward-backward charge asymmetry Ap B is defined from the lepton angular distribution with 
respect to the quark direction in the centre-of-mass frame as da/d cos 0* oc |(l + cos 2 9*) + Ap B cos9* . The 
lepton angle 9* in the dilepton center-of-mass frame can be calculated using the measured four momenta 
of the dilepton system; Ap B can then be determined with an unbinned maximum likelihood fit to the 
cos 9* distribution. 
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universe (> 1 Gyr). However, there are also tantalizing hints for the presence of an extra 
relativistic component, dubbed dark radiation. 

Taking these hints at face value, the most straightforward variation of standard ACDM 
is "extra" energy contributed by new relativistic particles U X." When the X's don't share 
in the energy released by e + e~ annihilation, it is convenient to account for the extra con- 
tribution to the SM energy density, by normalizing it to that of an "equivalent" neutrino 
species [101] 

7 

px = AN v p v = -AiV„p 7 , (49) 

o 

where p u is the energy density in neutrinos and p 7 is the energy density in photons (which 
by today have redshifted to become the CMB photons at a temperature of about 2.7 K). 
For each additional "neutrino-like" particle (i.e. any two-component fermion), if Tx = T UL , 
then AN U = 1; if X is a scalar (and Tx = T UL ), then AN V = 4/7. However, it may 
well be that the X's have decoupled even earlier in the evolution of the universe and have 
failed to profit from the heating when various other particle-antiparticle pairs annihilated (or 
unstable particles decayed). In this case, the contribution to AN^ from each such particle 
will be < 1 and < 4/7, respectively. The contribution of the 2.984 ± 0.009 neutrino species 
(measured from the width for invisible vv decays of the Z boson [102]) to N° s = N u + AN U 
is N v = 3.046; the small deviation from 3 is due to partial heating of neutrinos in the early 
universe by e + e~ annihilation, see e.g. [103]. 

The competition between gravitational potential and pressure gradients is responsible for 
the peaks and troughs in the CMB temperature angular power spectrum, see e.g. [104]. The 
redshift z eq of matter-radiation equality, 



1 + z, 



r 



8 111 



(50) 



affects the time (redshift) duration over which this competition occurs. Here, Q m h 2 is the 
total matter density (comprised, for nearly massless neutrinos, of baryons and CDM), h 
(H Q = lOO/i km/s/Mpc) is the normalized Hubble constant, and Vt^h 2 = 2.469 x 10 -5 is the 
present-day photon energy density. The primary effect of extra relativistic degrees of freedom 
on the CMB results essentially from changing the redshift of matter-radiation equality. If 
the radiation content is increased, matter-radiation equality is delayed, and occurs closer (in 
time and/or redshift) to the epoch of recombination. This implies the universe is younger 
at recombination with a correspondingly smaller sound horizon s*. Since the location of 
the n th peak in the angular power spectrum scales roughly as rmD^js* (where is the 
comoving angular diameter distance to recombination 11 ), if AN U > the peaks shift to 
smaller angular scales and with greater separation [104]. Therefore, the equality redshift is 
one of the fundamental observables that one can extract from WMAP data, mainly from 
the height of the third acoustic peak relative to the first peak. 



The angles on the sky are related to actual physical distance via the angular diameter distance d, defined 
as the ratio of the physical length (transverse to the line of sight) and the angle it covers d = A p h ys /$- 
Likewise, D = A c /$, where A c = (1 + z)A p h ys is the corresponding comoving length and z the redhsift; 
D = (l + z)d. 
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The variation in Nf reads [12] 

SN eS 5(tl h 2 ) 5z 

°J\ ~ 2.45 { " m 7 % > - 2.45 . (51) 

Nf n m h 2 i + z eq K J 

The latest distance measurements from the Baryon Acoustic Oscillations (BAO) in the dis- 
tribution of galaxies [16] and precise measurements of the Hubble constant Ho [17] provide 
an independent determination in the fractional error in Q m h 2 and allow a precise determi- 
nation of Nf. The parameter constraints from the combination of WMAP 7-year data, 
BAO, and H lead to Nf = 4.34±g;|| [12]. Similarly, a combination of BAO and H with 
data from the Atacama Cosmology Telescope (ACT) yields Nf = 4.6 ± 0.8 [105], whereas 
data collected with the South Pole Telescope (SPT) combined with BAO and Hq arrive at 
Nf = 3.86 ± 0.42 [106]. Although none of these measurements individually deviates from 
the standard value by more than about two standard deviations, they collectively rule out 
= 3.046 at the approximately 99% CL, and instead prefer roughly one extra effective 
neutrinos species [107]. 12 

The expansion rate of the universe at early times increases with the number of relativistic 
particle species in thermal equilibrium, and this in turn sets timescales for BBN [109- 
111]. One can then use the BBN yields of light nuclei to constrain the number of light 
species quantitatively. The nucleosynthesis chain begins with the formation of deuterium in 
the process p(n, 7)D. However, photo-dissociation by the high number density of photons 
delays production of deuterium (and other complex nuclei) until well after T drops below 
the binding energy of deuterium, Arj = 2.23 MeV. The number of photons per baryon above 
the deuterium photo-dissociation threshold, 7]~ 1 e~ AD ^ T , falls below unity at T ~ 0.1 MeV, 
where ?/ = %/n 7 ~ 5x 10 -10 is the baryon to photon number density. Nuclei can then begin 
to form without being immediately photo-dissociated again. Only 2-body reactions such as 
D(p, 7) 3 He, 3 He(D,p) 4 He, are important because the density is rather low at this time. 
Nearly all the surviving neutrons when nucleosynthesis begins end up bound in the most 
stable light element 4 He. Heavier nuclei do not form in any significant quantity both because 
of the absence of stable nuclei with mass number 5 or 8 (which impedes nucleosynthesis via 
n 4 He, p 4 He, or 4 He 4 He reactions) and the large Coulomb barriers for reactions such as 
T( 4 He, 7) 7 Li and 3 He( 4 He, 7) 7 Be. Hence the primordial mass fraction of 4 He, conventionally 
referred to as Y p , can be estimated by the simple counting argument 

>p = ^#. (52) 
p 1 + n/p v ' 

For T > 1 MeV, weak interactions were in thermal equilibrium, thus fixing the ratio of 
the neutron and proton number densities to be n/p = e~^^ T , where Q = 1.293 MeV 
is the neutron-proton mass difference. As the temperature dropped, the neutron-proton 
inter-conversion rate, T n ^ p ~ G^T 5 , fell faster than the Hubble expansion rate, H « 
a/ N(T) T 2 /M P1 (see e.g. [112]). Since N(T) counts the number of relativistic particle 
species determining the energy density in radiation, the neutron fraction n/p is directly 
sensitive to AN U . For standard ACDM, the freeze-out temperature of vl is 

/ 1 1/3 



FO 



Mpi G 



F 



1 MeV , (53) 



12 A more recent study seems to indicate 3.0 < Nf < 4.1 [108]. 
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yielding n/p ~ 1/6 and Y p ~ 0.25. 

The evidence for extra radiation from Y p data is, however, somewhat ambiguous. The 
observationally-inferred primordial fractions of baryonic mass in 4 He (Y p = 0.2472 ± 
0.0012 [113], Y p = 0.2477 ± 0.0029 [114], and Y p = 0.250 ± 0.004 [115]) have been con- 
stantly favoring Nf < 3 [116]. Unexpectedly, two recent independent studies deter- 
mined Y p = 0.2565 ± O.OOl(stat) ± 0.005(syst) [117] and Y p = 0.2561 ± 0.0108 [118]. For 
r n = 885.4 ± 0.9 s and r n = 878.5 ± 0.8 s, the updated effective number of light neutrino 
species is reported as Nf = 3.68±g;^ (2a) and Nf = 3.80±g;fg (2a), respectively [117]. 13 

As mentioned above, the primordial deuterium abundance depends not just on Nf but 
also on the cosmo logical baryon density, VL^h 2 . Prior to the precise inference of Qt>h 2 from 
CMB measurements, the strongest constraint on Nf came from using the deuterium-to- 
hydrogen number ratio D/H to restrict Vl^ti 2 and then exploiting the Nf dependence of the 
primordial helium mass fraction Y p . However, D /H has its own dependence on Nf [123, 124]; 
a strong external constraint on VL^h 2 allows BBN constraints on Nf that are independent 
of Y p . Since precise measurements of Y p are difficult, the constraint on Nf from D/H is 
found to be competitive with that from Y p . A recent analysis, which combines the CMB 
results with BBN theory and the observed D/H, suggests Nf = 3.90 ± 0.44 [125]. 14 

In summary, though uncertainties remain large, the most recent cosmological observa- 
tions show a consistent preference for additional relativistic degrees of freedom (r.d.o.f.) 
during BBN and the CMB epochs. We take these hints as motivation for the subsequent 
analysis, which consists of the following tasks: (1) to explain the dark radiation using the 
non-supersymmetric U(3)b x Sp(1)l x U(1)l X U(l)i R D-brane model, in which the addi- 
tional r.d.o.f. are the three flavors of light right-handed neutrinos which interact with the 
SM fermions via the exchange of heavy vector fields Z' and Z"; (2) to suppress the six addi- 
tional fermionic r.d.o.f. to levels in compliance with BBN and CMB. This is accomplished by 
imposing the decoupling of vrs from the plasma early enough so that they undergo incom- 
plete reheating during the quark-hadron transition; and late enough so as to leave an excess 
neutrino density suggested by the data. These requirements strongly constrain the masses 
of the heavy vector fields. Together with the couplings given in Table IV, the model is fully 
predictive, and can be confronted with dijet and dilepton data from LHC8 and, eventually, 
LHC14. 



B. Cosmology of Intersecting Branes 

The ensuing discussion will be framed in the context of a Z" which is mostly Jr, and we 
will comment on the case in which Z" is mostly B — L after presenting our results. 

We begin by first establishing, in a model independent manner, the range of decoupling 
temperatures implied by the BBN and CMB analyses. For the subsequent study, the physics 
of interest will be taking place at energies in the region of the quark-hadron transition, so 



For several years the Particle Data Group recommended r n = 885.7±0.8 s [119]. More recently, conflicting 
lifetimes r„ = 878.5 ±0.7 ±0.3 s [120] and r„ = 880.7 ± 1.3 ± 1.2 s [121] have been reported. The Particle 
Data Group now recommends a world average that includes the conflicting values, r„ = 881. 5± 1.5 s [122], 
with errors that have been inflated to reflect the discrepancy. 

The BBN calculations in this analysis includes updates of nuclear rates in light of recent experimental 
and theoretical information, with the most significant change occurring for the d(p, 7) 3 Hc cross section. 
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that we will restrict ourselves to the following fermionic fields, and their contribution to 
r.d.o.f.: [3u R ] + [3d R ] + [3s R ] + [3u L + e L + fi L } + [e R + fi R ] + [3u L + 3d L + 3s L ] + [3v R \. This 
amounts to 28 Weyl fields, translating to 56 fermionic r.d.o.f. 

Next, in line with our stated plan, we use the data estimate to calculate the range of 
decoupling temperature. The effective number of neutrino species contributing to r.d.o.f. 
can be written as N° s = 3[1 + (T UR /T UL ) 4 } ; therefore, taking into account the isentropic 
heating of the rest of the plasma between v R decoupling temperature T dec and the end of 
the reheating phase, 

(54) 

where T end is the temperature at the end of the reheating phase, and N(T) = r(T)(N B + ^NF) 
is the effective number of r.d.o.f. at temperature T, with JVb = 2 for each real vector field 
and Np = 2 for each spin-- Weyl field. The coefficient r{T) is unity for the lepton and 
photon contributions, and is the ratio s(T)/ssb for the quark-gluon plasma. Here s(T)(ssb) 
is the actual (ideal Stefan-Bolzmann) entropy. Hence N(T dec ) = 47.5 r(Td ec ) + 14.25. We 
take N(T cnd ) = 10.75 reflecting (e^ + e^ + e^ + e^z/ eL + z/ eR + ^ L + ^ R + z/ TL + z/ rR + 7 L + 7 R ). 
We consistently omit v R in considering the thermodynamics part of the discussion, but will 
include it when dealing with expansion. As stated in the introduction 

AM =/ °- 68 -o:35 M BBN + Y p 

u 10.90^ (la) CMB + BBN + H/D 1 1 

so the excess r.d.o.f. will lie within la of the central value of each set of observations if 
0.46 < AN U < 1.08. From Eqs. (54) and (55), the allowable range for N is 23 < N(T dcc ) < 
44. This is achieved for 0.18 < r(T dcc ) < 0.63. By comparing to Fig. 8 in Ref. [126], this can 
be translated into a temperature range 

175 MeV < T dcc < 250 MeV , (56) 

with the lower temperature coinciding with the region of most rapid rise of the entropy. Thus, 
the data implies that the v R decoupling takes place during the quark-hadron transition. 

We now turn to use our model in conjunction with the decoupling condition to constrain 
its parameters. To this end we calculate the interaction rate T(T) for a right-handed neutrino 
and determine T dec from the plasma via the prescription 

T(T dcc ) = H(T dcc ) . (57) 

Let f l L be a single species of Weyl fermion, representing the two r.d.o.f. {/£, f R }, where 
the superscript indicates bins i — 3,5. Similarly f R e {f l R , f l L }, for i — 1,2, 4, 6. Notice that 
the subscripts L, R denote the actual helicities of the massless particles in question, not the 
chirality of the fields. With this said, we may write the amplitude for f l L scattering 

m (u R { Pl )f L { P2 ) ^ Mp 3 )PM) = %HPBh^+iMpi)Mp*h^-i5)<P2)} ■ (58) 

The other 3 amplitudes are obtained by the crossing substitutions in the second square 
bracket; for scattering from 

fk -> v(p 2 ) 7^ (1 - 75) v(p 4 ) 

P R u(p 4 ) lfl (l + l5 )u(p 2 ) (59) 

fl "(P2) Ifi (1 + 7s) «(P4) ■ 
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The cross sections for the four scattering processes (no average over helicities) are 

o {vrA v R fi) = \<J [vrPr -+ y R f R ) = |^ (for bins i = 3, 5) (60) 

and 

o {ynfi -»> ^/i) = \°{vrSr -»> ^/«) = ^ — (forbinsz = 1,2,4,6) . (61) 

In addition to these scattering processes, the v R interacts with the plasma through the 
annihilation processes: vbPl — > fhfw f° r D i ns 2 = 3,5, and v>r&l — > PrPl-, f° r bins i = 
1, 2, 4, 6. These all yield cross sections 2G 2 s/ (37r) due to forward and backward suppression. 
Assuming all chemical potentials to be zero, the plasma will have an equal number density 
n(T) = 0.0913T 3 , for each fermion r.d.o.f. Thus, 

r scat (T) =n{T) l^a^VMMS , (62) 

where = 1 — cos6 l 12 is the Moller velocity, s = 2kik 2 (l — cos6 l 12 ) is the square of the 
center-of-mass energy, and Mi is the multiplicty of Weyl fields in each bin (e.g., for i = 
3, A3 = 3 + 2 = 5). The scattering cross section is given by 

4 2G 2 s 

af at = <r{y R f L -> v R f L ) + a{u R f R -> v R f R ) = -— *- for each i = 1, . . . , 6 ; (63) 

O 7T 



similarly, 



1 2G 2 i 

ar(s) = o{y R v L + = -— for each* = 1, . . . 6 . (64) 

O 7T 



Since s = 2k\k2(l — cos 9 12) and % = 1 — cos #12, we perform an approximate angular 
average ((1 — cos6 l i 2 ) 2 ) = 4/3, followed by a thermal averaging (2fci/c 2 ) = 2(3.15 2 T 2 ) to give 



-\scat 



(T) = (J^J ^2 (3.15T) 2 (0.0919T 3 ) (y^G 2 ^ ~ 2.05G 2 ff T 5 . (65) 



^From (63), (64), and (65), 

T ann (T) = Ir scat (T) ~ 0.50 G 2 ff T 5 . (66) 
Each of the Gi is given by the sum of the contributions from Z' and Z" exchange, 

C n' n' n" n" 

A—- — ^ 6 4. ^6 Hi /r>j\ 

The Hubble expansion parameter during this time is 

H(T) = 1.66 (A(T)) 1/2 T 2 /M P1 , (68) 



24 



Annihilation Only 
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FIG. 5: The shaded areas show the region allowed from decoupling requirements to accommodate 
CMB and BBN data. The hatched region indicates the masses excluded by the LHC8 dijet searches. 
The lower and upper shaded areas pertain to chemical and thermal equilibrium, respectively. These 
two estimates should serve to bracket the size of the actual effect. 

where Mp\ is the Planck mass. Since the quark-gluon energy density in the plasma has a 
similar T dependence to that of the entropy (see Fig. 7 in [126]), we take N(T) = 47.5 r(T) + 
19.5, so that H(T) = 10.3 T 2 /Mp\. (The first factor provides an average for r(T) over the 
temperature region, and we have now included the six vr r.d.o.f.) Since T oc T 5 and H ~ T 2 , 
it is clear that if at some temperature Td ec , H{Td ec ) = rj(T dec ), the ratio T/H will fall rapidly 
on further cooling. Thus from (57) and (68) the equation determining T dec depends on: (1) 
whether we need to preserve the absence of a chemical potential, or (2) whether we need 
simply to mantain physical equilibrium. The decoupling condition in these two cases is: (1) 
r ann (T dcc ) = H(T dec ) and (2) r scat (T dcc ) + r ann (T dcc ) = #(T dec ); or numerically: (1) 

0.50 G 2 ff T d 5 cc = 10.3 T d 2 cc /M pl =► T d 3 cc = 20.6 (G^Mp,)" 1 , (69) 

and (2) 

2.50 G 2 ff Tl c = 10.3 T 2 CC /M P1 =► T d 3 cc = 4.1 (G^Mn)" 1 . (70) 

T dcc as determined from these equations must lie in the band (56). 

Since all freedom of determining coupling constant and mixing angles has been exercised, 
there remains only constraints on the possible values of Mz> and Mz» ■ For high mass string 



25 



scales the contribution from Mz> to G e s is neglible. We find that for certain ranges of Mz" 
the decoupling of the v>rS occurs during the course of the quark-hadron transition, just so 
that they are only partially reheated compared to the u^s — the desired outcome. Since our 
aim is to match the data, which has lower and upper bounds on the neutrino "excess" , we 
obtain corresponding upper and lower bounds on the Z" gauge field mass. Roughly speaking, 
if decoupling requires a freezout of the annihilation channel (loss of chemical equilibrium), 
then 3.6 TeV < Mz» < 4.8 TeV. This range will be probed at LHC14. If thermal equilibrium 
via scattering is sufficient, then 5.4 TeV < Mz» < 7.4 TeV. 

Depending on the details of the string type model and M s some of the couplings may 
go up and some may go down, but the net result for G e s involving the product of all these 
couplings is virtually unchanged. Moreover, we have verified that if M s is pushed downwards 
to the TeV-scale region both Mz< and Mz" contribute to G e g and are within the LHC reach. 
A summary of LHC7 constraints and Mz< — Mz» mass regions consistent with CMB + BBN 
+ Yp + H/D data (within la) is encapsulated in Fig. 5. 

We comment briefly on the case in which Z" is mostly B — L. By comparing Tables IV 
and V it becomes evident that the Z" coupling to neutrinos is stronger when the extra gauge 
boson is almost diagonal in B — L. As a consequence, the allowed range of masses from 
decoupling requirements to accommodate CMB and BBN data is shifted to higher values: 
4.5 TeV < Mz" < 6.1 TeV if decoupling requires a freezout of the annihilation channel, and 
6.3 TeV < Mz" < 8.2 TeV if thermal equilibrium via scattering is sufficient. 

The first cosmology results from the Planck satellite anticipated in early 2013 would 
allow determination of iV° ff with a standard deviation of about 0.3, whereas the future 
Large Synoptic Survey Telescope (LSST) could determine N* G with a standard deviation of 
about 0.1 [127]. These observations when combined with future LHC results can directly 
test the viability of our model. 

V. SUPERSYMMETRIC EXTENSION 

When the string scale is at high energies, supersymmetry is in principle welcome for the 
hierarchy problem. Gauge bosons of the brane stacks belong then to J\f — 1 vector multiplets 
together with the corresponding gauginos, while at brane intersections chiral fermions belong 
to chiral multiplets denoted by their left-handed fermionic components Q, L, U c , D c , E c , N c , 
where the superscript c stands for the charged conjugate in the familiar notation. Moreover, 
in the (P, R) intersection, one should have the two usual Higgs doublets chiral multiplets 
Hi, H 2 with the quantum numbers of H* and H, respectively. Finally, the extra Higgs singlet 
H" becomes naturally the superpartner of the right-neutrino superfield N c . The Yukawa 
interactions (25) are now replaced by the superpotential: 

W Y = Y u QH 2 U c + Y d QH 1 D c + Y e LH 1 E c + Y N LH 2 N c . (71) 

On electroweak symmetry breaking, H 2 develops a VEV, as a result of which N c couples 
with vl to form a Dirac neutrino. Since superpotentials such as MN C N C or SN C N C are 
precluded by the U (l)x and U(l)i R gauge invariances, there seems no equivalent of the seesaw 
mechanism to generate the Weinberg term [128] which gives rise to Majorana neutrinos. 15 



However it is possible that D-branc instantons can generate Majorana masses for these perturbatively 
forbidden operators [129, 130]. 
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Here M is a Majorana mass matrix in flavor space and S is a gauge singlet. In addition, 
the existence of the VEV (N c ) breaks the U(1)l lepton gauge symmetry which allows the 
Z" to grow a mass. It also generates the i?-parity breaking term LH2, whose coefficient is 
subject to a variety of phenomenological constraints [131]. 

A superfield H" with Ir — L = +1 opposite to N c presents difficulties. A VEV for this 
version of H" serves equally well for the purpose of mass growth for Z" . However, its presence 
introduces a non-zero anomaly in B — L and Ir. The anomaly free status of Ir and B — L can 
be regained by introducing a fourth flavor N c . With this extension, the dimension 5 operator 
(N C H") 2 is permitted. This gives rise to a Majorana mass contrbibution oc v" /M s and to 
a pseudo-Dirac neutrino mass matrix [134, 135]. Present limits on pseudo-Dirac splittings 
arise from the solar and atmospheric neutrino measurements. Splitting of less than about 
10~ 12 eV 2 (for V\ and u 2 ) have no effect on the solar neutrino flux, while a pseudo-Dirac 
splitting of z/3 could be as large as 10~ 4 eV 2 before affecting the atmospheric neutrinos [136]. 
An even stronger bound emerges if we require the extra relativistic degrees of freedom not 
to exceed 1 as indicated by recent cosmological observations. To see this, we note that the 
effective thermalization of the right handed neutrinos can occur through mixing. This will 
occur if the oscillation length is less than horizon size during the CMB era. For a typical 
neutrino mass of ~ 0.1 eV, this requires that the Majorana mass is less than (9(10 -25 eV). 
At present we have no understanding of the origin of such a hierarchy (i.e. 10~ 13 beyond 
the ordinary suppression of the Yukawa), and as a consequence we discard the assignment 
Ir = L = +1 on phenomenological grounds. 

Like other broad frameworks for model-building, supersymmetric D-brane models do not 
lead uniquely to a single theory. 16 However, the conjectured models are rather rigidly con- 
strained, and lead to LHC predictions that are qualitatively different from the conventional 
minimal supersymmetric SM extensions [133]. 

We turn now to discuss some specifics of the SUSY extension to our analysis. The first 
and obvious change is the modification of the (3 functions for the running of the couplings. 
However, these changes will be minor: the phenomenological requirements at the TeV scale 
will effectively fix the U(l) couplings at that scale. Since unification is not a requirement 
of D-brane models, the coupling constants at the string scale will differ somewhat due to 
the change in the f3 functions, but string scale couplings do not alter our phenomenological 
predictions. The only caveat is to ensure that, as a result of the enhanced /3 functions, none 
of the couplings which comply with TeV data acquire non-pertubative components at the 
string scale. We have verified that the variation of the g[(M s ) parameter space is hardly 
noticeable. This gives scarcely any change in the production cross section and/or branching 
fractions, even in the extreme cases shown in Figs. 3 and 4, in which Z" is mostly diagonal 
in B — L or mostly diagonal in Ir. Furthermore, the milli-weak interactions required to 
explain the extra relativistic degrees of freedom during BBN and CMB epochs are largely 
independent of these changes. 

Much more serious considerations come to light in transcribing the low energy effective 
theory into a broken SUSY background. The technical problem arises most prominently in 
finding a broken SUSY framework that will accommodate the hierarchy between the mass of 
the Z and the mass of Z" . Breaking of the extra U(l) via the Biggs mechanism modeled on 
the radiative breaking of '577(2) xU(l) driven by a large top Yukawa coupling is not an option 



Some phenomenological aspects of the U(3)b x Sp(l)z x U(l)j R x U(1)l SUSY extension have been 
discussed in [132]. 
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in the present model. The introduction of an added D-term, a Fayet-Iliopoulos term, and an 
extended set of soft breaking masses, requires a sizable enlargement of the parameter space 
of the model. In order to incorporate this parameter space in a phenomenological study it 
is imperative to have additional experimental constraints on the SUSY spectrum. 

The approach we have taken here can be regarded as an effective theory with a new and 
novel phenomenology, as well as interesting theoretical characteristics (e.g., conservation of 
B to prevent proton decay and violation of L without Majorana masses). Of course such 
an effective theory requires a high level of fine tuning, which could be resolved in a more 
complete broken SUSY framework. However, we do not expect the phenomenology to differ 
in any substantial degree with the one presented in this paper. 

VI. CONCLUSIONS 

The main purpose of this paper has been to cast D-brane ideology in as bottoms-up, 
phenomenologically driven a way as possible. The energy scale associated with string physics 
is assumed to be near the Planck mass. To develop our program in the simplest way, we 
considered a minimal model with gauge-extended sector U(3)b x Sp(l)i x U(l)i R x U(1)l- 
The resulting U(l) content gauges the baryon number B, the lepton number L, and a 
third additional abelian charge Ir which acts as the third isospin component of an SU(2)r. 
Rotation of the U(l) gauge fields to a basis exactly diagonal in hypercharge Y and very 
nearly diagonal in (anomalous) B and (non-anomalous) Ir fixes all mixing angles and gauge 
couplings. The anomalous Z' gauge boson obtains a string scale Stiickelberg mass via a 4D 
version of the Green-Schwarz mechanism, TeV <C Mz> < M s < Mp\. To keep the realization 
of the Higgs mechanism minimal, we add an extra SU(2) singlet complex scalar, which 
acquires a VEV and gives a TeV-scale mass to the non-anomalous gauge boson Z". It is 
noteworthy that there are no dimension 4 operators involving H" that contribute to the 
Yukawa Lagrangian in our D-brane construct. This is very important since H" carries the 
quantum numbers of right-handed neutrino and its VEV breaks lepton number. However, 
this breaking can affect only higher- dimensional operators which are suppressed by the high 
string scale, and thus there is no phenomenological problem with experimental constraints 
for M s higher than ~ 10 14 GeV. Since all freedom of determining coupling constant and 
mixing angles has been exercised, there remains only constraints on the possible value of 
Adz" ■ We have shown that Mz» ~ 3 — 4 TeV saturates current limits from the CMS and 
ATLAS collaborations. We have also shown that for Mz» < 5 TeV, LHC14 will reach 
discovery sensitivity > 5a. 

Armed with our D-brane construct, we developed a dynamic explanation of recent hints 
that the relativistic component of the energy during the CMB and BBN epochs is equivalent 
to about 1 extra Weyl neutrino. Requiring that the B — L current be anomaly free implies 
existence of 3 right-handed Weyl neutrinos. The task then reverts to explain why there 
are not 3 additional r.d.o.f. We showed that for certain ranges of Mz» the decoupling of 
the vrS occurs during the course of the quark-hadron crossover transition, just so that 
they are only partially reheated compared to the z^'s — the desired outcome. Roughly 
speaking, if decoupling requires a freezout of the annihilation channel (loss of chemical 
equilibrium), then for a Z" which is mostly Ir, 3.6 TeV < Mz» < 4.8 TeV, whereas for 
a Z" which is mostly B — L, 4.5 TeV < Mz» < 6.1 TeV. This range will be probed at 
LHC14. If thermal equilibrium via scattering is sufficient, for a Z" which is mostly Ir, 
5.4 TeV < M z » < 7.4 TeV, and for a Z" which is mostly B - L, 6.3 TeV < M z » < 
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8.2 TeV. To carry out this program, we needed to make use of some high statistics lattice 
simulations of a QCD plasma in the hot phase, especially the behavior of the entropy during 
the confinement- deconfinement changeover. Interestingly, the behavior of the trace anomaly 
(shown in Fig. 15 of [126]), which is very sensitive to the nature of the crossover region, 
shows a sharp peak at 200 MeV and our range for T^ ec straddles this region. 

Throughout this paper we remained agnostic with respect to SUSY breaking and the 
details of the low energy effective potential. However, we do subject the choice of quantun 
numbers for H" to the stringent holonomic constraints of the superpotential at the string 
scale. This forbids the simultaneous presence of scalar fields and their complex conjugate. 
As an illustration, if the quantum numbers of H" are those of N^, then higher dimensional 
operators such as NrN^H" 2 , which can potentially generate a Majorana mass, are absent. 
Because of holonomy this absence cannot be circunvented by including N rNrH"* 2 . 

In summary, we have studied the U(l) phenomenology of D-brane models endowed with 
a high mass string scale. We have incorporated some elements of SUSY, discussing evo- 
lution of the gauge couplings to the string scale and enforcing the holonomic constraints 
on the superpotential. We have shown that LHC8 data set upper limits on the mass of 
the Z" gauge boson: Mz" ^3 — 4 TeV. We have also shown that Z" milli-weak inter- 
actions, which are within reach of LHC14, could play an important role in observational 
cosmology. It is important to stress that the Z" production cross section and its branching 
fractions are universal and have been evaluated in a parameter-free manner. Therefore, the 
U(l) phenomenology presented in this paper is completely independent of the details of 
the compactification scheme, such as the configuration of branes, the geometry of the extra 
dimensions, and whether the low energy theory is supersymmetric or not. 
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